Chaos in resonant-tunneling superlattices 
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Spatio-temporal chaos is predicted to occur in n-doped semiconductor superlattices with sequential 
resonant tunneling as their main charge transport mechanism. Under dc voltage bias, undamped 
time-dependent oscillations of the current (due to the motion and recycling of electric field domain 
walls) have been observed in recent experiments. Chaos is the result of forcing this natural oscillation 
by means of an appropriate external microwave signal. 
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Nonlinear oscillations and chaos have been predicted 
in systems where a predominantly quantum dynamics is 
corrected by mean-field nonlinear terms due to collective 
interactions (Hartree)lil or to interactions with classical 
subsystems having a widely different time scako These 
phenomena are different from quantum chaos,Q i.e., the 
behavior of quantum systems whose classical counterpart 
is chaotic. So far chaotic oscillations have been predicted 
for systems with few degrees of freedom and experimen- 
tal evidence is scarce. In this paper we predict chaotic 
behavior with loss of spatial coherence in a system with 
many degrees of freedom for which the main transport 
mechanism is resonant tunneling: a weakly-coupled mul- 
tiquantum well superlattice (SL). As far as we know these 
are the first results on chaotic behavior in SLs. In con- 
trast to unbiased triple-well heterostructure,El our SL is 
subject to external dc+ac bias, as it was the case with 
the two-level system of Ref. 2. 

Very recently time-dependent oscillations of the cur- 
rent on GaAs/AlAs SLs subject to dc voltage bias have 
been found.Q The oscillations are damped for undoped 
photoexcited SLs and undamped for doped SLs without 
photoexcitation.El For large values of the photoexcitation 
or the doping, there is stable formation of stationary 
electric field domains leading to the well-known oscilla- 
tory I-V characteristic.ErEl'El According to a discrete drift 
modeljEl the current oscillations are caused by the cre- 
ation, motion and recycling of domain walls separating 
two electric field domains.EHI This situation is reminiscent 
of that found in bulk semiconductor devices with nega- 
tive differential resistance (NDR), where dc voltage bias 
gives rise to high-field domain dynamics and the well- 
known Gunn oscillations Oll3 A significant difference is 
that the space charge waves are dipoles in the Gunn os- 
cillations and charge monopolcs in the SL current oscil- 
lations. Another is that the Gunn waves are generated 
close to the injecting contacj—whereas the domain walls 
appear clearly inside the SL.E3 

Having found a system with a natural oscillation due to 
traveling-wave motion, it is natural to ask whether har- 
monic forcing would lead to chaos with spatial structure. 
The answer is affirmative. This is also the case for the 
periodically driven Gunn diode studied by Mosekilde et 



perimental studies have been carried out by Kahn 
et alE3 on ultrapure p-Ge, where the NDR is caused by 
negative differential impurity impact ionization,E_3 and 
the transition to a chaotic attractor with loss of spatial 
coherence has been observed. 

We consider a set of weakly interacting quantum wells 
(QWs) characterized by average values of the electric 
field £i(t), and the electron density rij(t), with i = 
1,...,N denoting the QW index. This mean-field-like 
approach is often justified because the relevant time scale 
for the oscillations (~ 0.1 /is) is much larger than those 
for the tunneling process between adjacent QWs (~ 1 ns) 
and the relaxation from excited levels to the ground state 
within each QW (~ 1 ps).B The one-dimensional equa- 
tions governing the dynamics of the system are the Pois- 
son equation averaged over one SL period I, Ampere's 
equation for the balance of current density and the volt- 
age bias condition: 

y (St- £i-x)= & - (m-N D ) (1) 

e ^ + eriiV^i) = J (2) 
at 

N 

lY J ^ = V(t). (3) 

i=l 

Here e, e, and No are the average permittivity, 
the electron— charge, and the average doping density, 
respectively.E3 The total current density J(t) is the sum 
of the displacement current and the electron flux due 
to sequential resonant tunneling eriiv(£i). The effective 
electron velocity v{£) exhibits maxima at the resonant 
fields for which the adjacent levels of neighboring QWs 
are aligned Jill as shown in the inset of Fig. [I]. The volt- 
age V(t) in (|^) is the sum of a dc voltage V& and an ac 
microwave signal of relative amplitude A and driving fre- 
quency f d : V(t) = V b {l + Asm(2irf d t)}. The boundary 
condition e(£±—£o)/(el) = ni — No = S allows for a small 
negative charge accumulation in the first well, which is 
taken to be 5 ~ 10 -3 x Nr>- The physical origin of 8 is 
that the n-doped SL is typically sandwiched between two 
n-doped layers with as excess of electrons, thereby form- 
ing a n + -n-n + diode.u Then some charge will be trans- 
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ferred from the contact to the first QW creating a small 
dipole field that will cancel the electron flow caused by 
the different concentration of electrons at each side of the 
first barrier. 

To study our equations, it is convenient to render them 
dimensionless, using the characteristic physical quanti- 
ties. As the unit of the electric field E = £/£x-2 we 
adopt the typical electric field strength needed to align 
the first and the second electron subbands of the nearest 
QWs £1-2 ■ The dimensionless velocity v(E) is obtained 
normalizing v(£) by its value at £ = £1-21 where it has a 
local maximum due to resonance in tunneling (sec inset of 
Fig. [I). The others dimensionless quantities are defined 
as follows: the donor concentration v = el No / e£x-2\ 
the time r = i/ttun where itun = l/v(£\-2) is the char- 
acteristic tunneling time; the dc bias V = V\,/£x-2lN; 
the ac bias amplitude a = AV; the driving frequency 
u> = 2'Kfdtt V Ln- By time-differentiating ([}]) and using 
Ampere's law (||), one can express the current density 
as 



e dV 

iN~dT 



e 
N 



N 



(4) 



and after its substitution into (^|) we obtain a system of 
N equations for the electric field profiles: 
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— y(Ei) [Ei — Ei-i + v\ + acjcos(wr), (5) 

which we solved numerically by the fourth-order Runge- 
Kutta method with the boundary condition Eq = Ex — vb 
and initial conditions Ei(0) = V, Vi. As an example, 
we consider a GaAs/AlAs SL at T=5 K with ^=40, 
1=13 nm, £i_ 2 «10 5 V/cm, N D sal .1 5x f 17 c m - 3 , for 
which undamped time-dependent oscillations of the cur- 
rent were first observed.El One get v w 0.1, t tun ~ 2.7 ns, 
and we take V=1.2 (corresponding to Vb ~7.8 V). The 
main features of our numerical results are as follows. 

For the dc case (a=0) undamped time-periodic current 
oscillations of frequency /o ftj-S MHz, (in excellent agree- 
ment with the observed valuetj) set in after a transient pe- 
riod. The electric field and charge profiles corresponding 
to one period of the current oscillation are similar to those 
found in the undoped case (cf. Figure 9 in Ref. 9). During 
each period a domain wall (charge monopole) is formed 
inside the SL. It then moves towards the corresponding 
contact. Depending on the applied voltage, it may or may 
not reach the end of the SL before it dissolves and a new 
monopole is formed starting a new period of the oscilla- 
tion. Since in a 40-period SL the wall moves only a few 
QWs, it may be hard to distinguish this monopole recy- 
cling from a true spatial oscillation of a single monopole. 
Simulations of longer SLs (N >100), however, clearly 
show monopole recycling with two monopoles coexisting 
during some part of one current oscillation period. The 



frequency of the oscillation is mainly determined by the 
number of QWs the monopole moves across (which in- 
creases with N) and by the average drift velocityJl3 

Now consider the ac case. We start with a uniform 
initial field profile and solve the equations for dc bias. 
After a short transient, the self-sustained oscillations set 
in and we switch on the ac part of the bias. Our main 
result is that the competition between the natural oscil- 
lation due to monopole dynamics and the forcing gives 
rise to narrow windows of spatio-temporal chaos for ap- 
propriate values of V, a and u>, and that the size and rich- 
ness of these windows increases with N. Let us fix the 
ratio between the natural frequency fo and the driving 
frequency fd at the golden mean (a/5 - l)/2. To detect 
and visualize the chaotic regions in parameter space, we 
need to define a Poincare mapping. The current is a good 
measure of the amplitude (norm) of the solutions which is 
illustrated by the use of current vs voltage characteristics 
as bifurcation diagrams £3 Denoting the period of the ac 
bias Td = 2tt/uj, we adopt as Poincare mapping (for each 
value of a) the current at times r m = mTd, m = 0, 1, . . . 
(after waiting enough time for the transients to have de- 
cayed). The result is the bifurcation diagram in Figj^. 
Notice the period-doubling sequences that point to the 
existence of chaos near their accumulation points. There 
we have computed the largest Lyapunov exponent and 
found it positive, which confirms chaos within the win- 
dows marked by arrows in FigJ^. The chaotic regions 
are interspersed with locking to periodic and quasiperi- 
odic regimes. Quasiperiodic routes to chaos have been 
found at the first and last windows marked with arrows 
in Figjl} Notice that the period-2 orbits span the widest 
parameter region from the narrow chaotic band around 
a ps 0.01 up to the next chaotic region at a ss 0.09. For 
a > 0.145 the solution is attracted to the period-1 orbit 
with the driving frequency fd- 

More insight into the transition between chaotic and 
non-chaotic regions of the bifurcation diagram can be 
obtained by sweeping-down through it, as follows. We 
set ao=0.16 (the leftmost value in Fig. [j| where there is 
pcriod-1 locking) and take an initial field profile {Ei} = 
V. We now solve the problem (||), compute J m = J (mTd) 
and find out when the solution is periodic within a 10 -5 
accuracy. At that time, we stop the simulation, depict 
all J m correspondent to the period, store the resulting 
electric field profile, and use it as initial condition for an- 
other simulation with a = ao — Aa, Aa = 2 x 10~ 4 . By 
repeating this process, we arrive at chaotic or quasiperi- 
odic regions of a, where the integration is stopped at a 
much larger time mTd, with m=2000, and only the last 
1000 points are depicted, thereby eliminating transients. 
Fig.|l| is the result of such a sweep-down run. The tran- 
sition points between different regions in the bifurcation 
diagram are found to be different when a sweep-up run 
is made, demonstrating thus hysteresis. 

To illustrate the spatially chaotic nature of the solu- 
tions we pick arbitrarily two far-away QWs and depict 
the simultaneous values of the electric field at them af- 
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ter each period of the driving force The resultant 
attractors for several amplitudes a are presented in Fig. 
pi The first example [Fig. 0(a)] shows a chaotic attrac- 
tor with layered structure and variation in the density 
of points. The chaotic attractor of Fig. ||(b) has several 
separate branches almost continuously filled. Finally Fig. 
|](c) corresponds to a value of a on the quasiperiodic re- 
gion. The closed loop with periodic pattern [see Fig. 
||(d)] indicates quasiperiodicity: the orbit fills the attrac- 
tor (torus), never closing on itself. The points in the 
chaotic attractors cluster with varying density on differ- 
ent regions, which means that they can be characterized 
by their multifractal dimension D q (calculated by the 
method of Ref. 18). We have found that for the partic- 
ular case of Fig. ||(a), D q decreases from D-oo « 1.57 
to D +OQ sa 0.72, going through the capacity dimension 
D « 1.18. 

Loss of spatial coherence in the chaotic regime is easily 
found in long SLs. Density plots for the electron con- 
centration (dark regions mean high electron density) in 
a 200-pcriod SL with natural frequency ~ 1.9 MHz are 
shown in Fig. ^. The density plots show the transition 
from periodic [Fig. |^(a)] to chaotic domain wall dynamics 
with loss of spatial coherence [Fig. ||(b)]. Under dc + ac 
driving with golden mean frequency ratio, nucleation of 
monopole wavefronts occurs more frequently: in addition 
to long- living waves travelling over almost the whole SL, 
there are short-living waves. The two types of waves are 
distributed chaotically in space and nucleated both at the 
beginning and deep inside the SL. Detailed calculations 
show coexistence of up to three monopoles connecting 
four electric field domains during certain time intervals. 
We thus have that the loss of spatial coherence is due to 
chaotic domain- wall dynamics, as seen in Fig. ||(b). 

In conclusion, spatio-temporal chaos is expected to oc- 
cur in weakly-coupled semiconductor superlattices (with 
sequential resonant tunneling as the main transport 
mechanism) under appropriate dc+ac voltage bias. This 
prediction should be experimentally testable in currently 
available n-doped GaAs/AlAs samples forming n + -n-n + 
diodes .B 
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FIG. 1. Bifurcation diagram of the current at times mTd 
(for sufficiently large m) versus the driving-force amplitude, 
for a 40-period SL and the golden-mean ratio between natural 
and driving frequencies. Windows of chaotic solutions are 
marked by arrows. Inset: Dimensionless velocity as a function 
of the electric field. The point indicates the electric-field value 
corresponding to the dc bias V = 1.2 used in the calculations. 
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FIG. 2. Electric field values E$ and (at two different 
QWs of the 40-period SL) measured at times mTd for: (a) 
a=0.088; (b) a=0.101; (c) a=0.14. (d) is a blow-up of the 
small region inside the rectangle in (c). 10 4 successive points 
have been used to depict each attractor. 
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FIG. 3. Density plots for the electron concentration in a 
200-period SL for: (a) pure dc bias; and (b) dc + ac bias 
with a=0.112. x- and y-axes correspond to time (in units of 
16 ttun) and QW index respectively. Darker regions indicate 
high electron density and mark the location of the domain 
boundary for each time. 
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